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Systems of interacting quantum spins show a rich spectrum of quantum phases and display in¬ 
teresting many-body dynamics. Computing characteristics of even small systems on conventional 
computers poses significant challenges. A quantum simulator has the potential to outperform stan¬ 
dard computers in calculating the evolution of complex quantum systems. Here, we perform a digital 
quantum simulation of the paradigmatic Heisenberg and Ising interacting spin models using a two 
transmon-qubit circuit quantum electrodynamics setup. We make use of the exchange interaction 
naturally present in the simulator to construct a digital decomposition of the model-specific evolu¬ 
tion and extract its full dynamics. This approach is universal and efficient, employing only resources 
which are polynomial in the number of spins and indicates a path towards the controlled simulation 
of general spin dynamics in superconducting qubit platforms. 


Quantum simulations using well controllable quan¬ 
tum systems to simulate the properties of another less 
tractable one [1, 2] are expected to be able to predict 
the properties and dynamics of diverse systems in con¬ 
densed matter [3, 4], quantum chemistry [5] and high 
energy physics [6, 7]. In particular, quantum simulations 
are expected to provide new insights into open problems 
such as modeling high-Tc superconductivity [8], ther- 
malization [9] and non-equilibrium dynamics [10]. Up 
to now, several prototypical quantum simulations have 
been proposed and realized in trapped ions [11], cold 
atoms [12], and quantum photonics [13]. Examples in¬ 
clude spin models [14-16], many-body physics [17], and 
relativistic quantum mechanics [18]. In the field of super¬ 
conducting circuits quantum simulations are still in their 
infancy [19]. Topological properties [20, 21] have been 
simulated recently, as have been fermionic models [22]. 

Quantum simulators are typically classified into two 
main categories, namely, analog and digital ones. Ana¬ 
log quantum simulators are designed to display intrin¬ 
sic dynamics which are equivalent to those of the sim¬ 
ulated system. While this approach is not universal it 
features control of the relevant Hamiltonian parameters 
better than in the system to be simulated. Instead, dig¬ 
ital quantum simulators [2] can reproduce the dynamics 
of a quantum system via a universal digital decompo¬ 
sition of its Hamiltonian H = Hk into efficient el¬ 
ementary gates realizing H^. This approach is based 
on the Suzuki-Lie-Trotter expansion of the time evolu¬ 
tion U(t) = e~ lHt = lim n _ ) . 00 (n^ =1 e - lH kt/n^n anc [ was 
recently demonstrated experimentally in a trapped-ion 
digital quantum simulator [15]. 

Here we demonstrate digital quantum simulation of 
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spin systems [16] in an architecture known as circuit 
quantum electrodynamics (QED) [23]. 

Our experiments are carried out with two supercon¬ 
ducting transmon qubits [24] coupled dispersively to a 
common mode of a coplanar waveguide resonator (see 
Appendix A for the device layout and setup diagram). 
We operate the circuit at 30 mK in a dilution refriger¬ 
ator. The qubits Q1 and Q2 interact with a coplanar 
waveguide resonator with a fundamental resonance fre¬ 
quency at 7.14 GHz which serves both as a quantum 
bus [25] and for readout [26]. 

The natural two-qubit interaction is the XY exchange 
coupling [25] H* y 2 = ^ (erferf + cr y a y ) mediated by vir¬ 
tual photons in a common cavity mode, which we also 
refer to as the XY interaction. Here, a^ ,y are the Pauli 
operators acting on qubit i and J denotes the effective 
qubit-qubit coupling strength [27]. The XY interaction 
is activated by tuning the transition frequency of qubit 
Q1 (5.44 GHz) into resonance with qubit Q2 (5.24 GHz) 
for a time r using nanosecond time scale magnetic flux 
bias pulses [28] (see Appendix B). When the qubit tran¬ 
sition frequencies are degenerate, the resonator-mediated 
coupling strength is spectroscopically determined to be 
J = —40.4 MHz. To make the presentation of the sim¬ 
ulation results independent of the actual J, we express 
the interaction time r for a given J in terms of the ac¬ 
quired quantum phase angle 2\J\r. In our setup, the 
action of the XY gate (Fig. la) is characterized by full 
process tomography for a complete set of 16 initial two- 
qubit states and a series of 25 different interaction times 
r finding process fidelities no lower than 89% (see Ap¬ 
pendix D). 

In Fig. 2a,b we present non-stationary spin dynamics 
under the XY exchange interaction for a characteristic 
initial two-qubit state |t)(|t) + \\))/V2 with spins point¬ 
ing in perpendicular directions along +z and +x, respec¬ 
tively. During the XY interaction, the state of one spin is 
gradually swapped to the other spin and vice versa with 



2 


(a) 



(b) 


XY 


XZ 


YZ 




FIG. 1 . (a) Circuit diagram to characterize the XY exchange 
interaction on the qubits Q1 and Q2 symbolized by the ver¬ 
tical line (x) which is activated for a time r. To perform 
standard process tomography of this interaction, separable 
initial states are prepared using single-qubit rotations R^^ 
(green) in the beginning and the final state is characterized 
using single-qubit basis rotations R\°™ and joint two-qubit 
readout (yellow), (b) Digital quantum simulation of the two- 
spin Heisenberg (XYZ) interaction for time r. The first step 
after state-preparation is to apply the XY gate for a time r 
(dashed box labeled as XY). In the second and third steps 
(dashed boxes with labels XZ and YZ), XZ and YZ gates are 
realized using single-qubit rotations Rt^j 2 (blue) by an angle 
±7r/2 about the x or y axis transforming the basis in which 
the XY gate acts, (c) Protocol to decompose and simulate 
Ising spin dynamics in a homogeneous transverse magnetic 
field. The circuit between the bold vertical bars with two 
dots is repeated n times, invoking each XY and phase gates 
for a time r/n. See text for details. The actual pulse scheme 
is provided in Appendix C. 


a phase angle of n/2. This corresponds to the iSWAP 
gate [29]. As a consequence, the measured Bloch vec¬ 
tors move along the YZ and XZ planes. For a quantum 
phase angle of 2| J\r = i r they point along the +y and 
+z directions respectively in good agreement with the 
ideal unitary time evolution indicated by dashed lines in 
Fig. 2a,b. We also find that the two-qubit entanglement 
characterized by the measured negativity [30] of 0.246 
is close to the maximum expected value of 0.25 for this 
initial state at a quantum phase angle of 7t/2. As a con¬ 
sequence the Bloch vectors do not remain on the surface 
of the Bloch sphere but rather lie within the sphere. 

The anisotropic Heisenberg model describes spins in¬ 
teracting in three spatial dimensions 

H xyz = ^2(J x a?(Tj + Jy<J V iCF V j + ( 1 ) 

(i,j) 

where the sum is taken over pairs of neighbouring spins i 
and j. J x , J y and J z are the couplings of the spins along 
the x, y and z coordinates, respectively. Since it does not 


occur naturally in circuit QED we decompose the Heisen¬ 
berg interaction into a sequence of XY and single-qubit 
gates, as shown in Fig. lb. We combine three succes¬ 
sive effective XY, XZ and YZ gates derived from the XY 
gate by basis transformations [16] to realize the isotropic 
Heisenberg model with J x = J y = J z = J versus in¬ 
teraction time r. Since the XY, XZ and YZ operators 
commute for two spins the Trotter formula is exact after 
a single step. 

To compare the Heisenberg (XYZ) interaction with the 
XY exchange interaction we have prepared the same ini¬ 
tial state as presented in Fig. 2a,b. The isotropic Heisen¬ 
berg interaction described by the scalar product between 
two vectorial spin 1/2 operators preserves the angle be¬ 
tween the two spins. As a result, the initially perpendicu¬ 
lar Bloch vectors of qubits Q1 and Q2 remain perpendic¬ 
ular during the interaction (Fig. 2c) and rotate clockwise 
along an elliptical path that spans a plane perpendicu¬ 
lar to the diagonal at half angle between the two Bloch 
vectors (Fig. 2c). 

In accordance with theory, the XYZ interaction leads 
to a full SWAP operation for a quantum phase angle 
of 21 J\t = 7r/2 where the Bloch vectors point along the 
+x and +z directions. For the given initial state, we 
observed a maximum negativity of 0.210 close to the ex¬ 
pected value of 0.25 for the Heisenberg interaction at 
a quantum phase angle of 2| J\r = 7r/4. As for the XY 
interaction we have characterized the Heisenberg interac¬ 
tion with standard process tomography finding fidelities 
above 82% for all quantum phase angles 2\J\r. 

Next, we consider the quantum simulation of the Ising 
model with a transverse homogeneous magnetic field 

+ (2) 

(bi) i 

where the magnetic field B pointing along the z axis is 
perpendicular to the interaction given by JcrfcrJ. Since 
the two-spin evolution (Fig. lc) is decomposed into two- 
qubit XY and single-qubit Z gates which do not com¬ 
mute, the transverse field Ising dynamics is only recov¬ 
ered using the Trotter expansion in the limit of a large 
number of steps n for an interaction time of r/n in each 
step. To realize the Ising interaction term using the ex¬ 
change interaction, the XY gate is applied twice for a 
time r/n, once enclosed by a pair of 7r pulses on qubit 
Ql. This leads to a change of sign of the term which 

thus gets canceled when added to the bare XY gate. The 
external magnetic field part of the Hamiltonian is re¬ 
alized as single-qubit phase gates Rf which rotate the 
Bloch vector about the £ axis by an angle </> = Br/n per 
Trotter step. These gates are realized by detuning the 
respective qubit by an amount S from its idle frequency 
corresponding to an effective 5-field strength of B = 2i tS. 

We experimentally simulate the non-stationary dy¬ 
namics of two spins in this model for the initial state 
It)(It) — i\\))/V2 which is well-suited to assess the simu¬ 
lation performance. In Fig. 3a expectation values for the 
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FIG. 2. (a) Experimentally determined coordinates of the Bloch vectors during exchange (XY) interaction represented by small 
red (Ql) and blue (Q2) points are compared to the ideal paths shown as dashed lines in the XY model. The ideal paths are 
in the YZ and XZ planes shown as blue and red planes intersecting the Bloch sphere. The time evolution is indicated by 
the saturation of the colors as the quantum phase angle 2| J\r advances from 0 (saturated) to 7r (unsaturated), (b) Measured 
expectation values of the Pauli operators oyfor the qubits Ql (red points) and Q2 (blue points), respectively, for the XY 
interaction as a function of the quantum phase angle 2| J\r along with the ideal evolution (dashed line), (c) Evolution of the 
Bloch vector for the quantum simulation of the isotropic Heisenberg interaction vs. quantum phase angles from 0 to 37r/4. The 
path of the Bloch vectors of the qubits Ql and Q2 spans the plane indicated by the rectangular sheets intersecting the Bloch 
spheres, (d), As in panel b for the Heisenberg interaction. 


digital simulation of the erf 2 -components of the two spins 
are shown, as well as the two-point correlation function 
(erferf). The erf 2 -components of the spins represented 
by the red and blue datasets in Fig. 3a, respectively, os¬ 
cillate with a dominant frequency component of 2J due 
to the presence of the interaction term oc erf erf. Like¬ 
wise, the XX correlation (erferf) represented by the yel¬ 
low dataset in Fig. 3a is non-stationary and oscillates at 
rate 2 y/B 2 + J 2 = 2^/lOJ « 6.3 J due to the presence of 
a magnetic field of strength B = 3 J. The evolution of the 
measured final state shows agreement with a theoretical 
model (solid lines in Fig. 3a) which takes into account 
dissipation and decoherence with deviations being dom¬ 
inated by systematic gate errors (see Appendix E). 

In Fig. 3b the fidelity of the simulated state is com¬ 
pared to the expected state at characteristic quantum 
phase angles both for the experimental realization (col¬ 
ored bars) and the ideal Trotter approximation (wire 
frames) after the nth step. In an ideal digital quantum 
simulator the theoretical fidelity (wire frame) converges 
for an increasing number of steps n (Fig. 3b). The ex¬ 
perimental fidelity, however, reaches a maximum for a 
finite number of steps (Fig. 3b) after which it starts to 


decrease due to gate errors and decoherence [16]. As ex¬ 
pected, the Trotter approximation converges faster for 
smaller quantum phase angles 2|J|r. For 2\J\r = 7r/4 
the peak experimental fidelity (Fig. 3b) of 98.3 % is al¬ 
ready observed for n = 1, whereas for 2| J\r = 3 tt/2 the 
optimum of 80.7% is observed for n = 5. 

In future experiments, transmission line resonators 
may provide a means to design multi-qubit devices with 
non-local qubit-qubit couplings that directly reflect the 
lattice topology of spin systems such as frustrated mag¬ 
nets. Moreover, the incorporation of cavity modes as 
explicit degrees of freedom in the simulated models [31], 
following an analog-digital approach, and the integration 
of optimal control concepts, will be instrumental to scale 
the system to larger Hilbert-space dimensions. With this, 
the circuit QED architecture offers considerable poten¬ 
tial for surpassing the limitations of classical simulations, 
which can be facilitated by using efficient digital decom¬ 
positions of spin Hamiltonians as pursued in this work. 
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FIG. 3. (a) Digital quantum simulation of the Ising model with transverse homogeneous magnetic field using 1 to 3 Trotter 
steps. Shown are the ^-components (erf) of qubit Q1 (red) and (erf) of qubit Q2 (blue) and the two-point correlation function 
in the ^-direction (erferf) (yellow points) of the spins as a function of the quantum phase angle 2\J\r for the initial state 
It)(It) — i\i))/V2 and a magnetic field strength B = 3J. Theoretically expected results take systematic phase offsets and finite 
coherence of the qubits into account (solid curves). The ideal dynamics are obtained from the time-dependent Schrodinger 
equation for the Ising Hamiltonian (dashed lines), (b) Fidelity with respect to the exactly solved Ising model for displayed 
quantum phase angles of the final state after ideal unitary evolution in the simulation protocol for n Trotter steps (wire frames) 
and experimentally obtained final state (colored bars). 
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Appendix A: Chip architecture and measurement 
setup 

The present experiment was performed using two su¬ 
perconducting transmon [24] qubits Q1 and Q2 and 
one coplanar waveguide resonator R1 on a microchip 
(Fig. Al). The resonator R1 has a fundamental reso¬ 
nance frequency of v r = 7.14 GHz. From spectroscopic 
measurements we have determined the maximum transi¬ 
tion frequencies ^ max = {5.55, 5.53} GHz and charging 
energies Ec/h ~ {260, 260} MHz of the qubits Q1 and 
Q2, respectively, where h is the Planck constant. The 
qubits Q1 and Q2 are coupled to resonator R1 with cou¬ 
pling strengths g/2n ~ {120, 120} MHz. For this exper¬ 
iment the qubit transition frequencies in their idle state 
were offset to v = {5.440, 5.240} GHz by applying a con¬ 
stant magnetic flux threading their SQUID loops with 


miniature superconducting coils mounted underneath the 
chip. At these idle frequencies, the measured energy re¬ 
laxation and coherence times were T\ = {7.1, 6.7} fi s and 
X 2 = {5.4, 4.9} fi s, respectively. The transition frequen¬ 
cies of the qubits Q3 and Q4 were tuned to 4.5 GHz and 
6.1 GHz such that they do not interact with Q1 and Q2 
during the experiment. 

A schematic diagram of the measurement setup is 
shown in Fig. A2a. To realize two-qubit XY gates and 
single-qubit phase gates (Z), controlled voltage pulses 
generated by an arbitrary waveform generator (AWG) 
are used to tune the flux threading the SQUID loop of 
each qubit individually using flux bias lines [28]. The 
single-qubit microwave pulses (X,Y) are generated us¬ 
ing sideband modulation of an up-conversion in-phase 
quadrature (IQ) mixer (Fig. A2b) driven by a local 
oscillator (LO) and modulated by an arbitrary wave¬ 
form generator (AWG). The same up-conversion LO is 
used for the microwave pulses on both qubits to mini¬ 
mize the phase error introduced by phase drifts of mi¬ 
crowave generators. We have used a quantum-limited 
parametric amplifier (PA) to amplify readout pulses at 
the output of R1 (Fig. A2c). Here the Josephson junc¬ 
tion based amplifier in form of a Josephson parametric 
dimer (JPD) [32] is pumped by a strong pump drive 
through a directional coupler (D). To cancel the pump 
leakage, a phase (0) and amplitude (A) controlled mi¬ 
crowave cancelation tone is coupled to the other port 
of the directional coupler (D). Three circulators (Cl-3) 
were used to isolate the sample from the pump tone. A 
circulator (C4) at base temperature followed by a cav- 
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FIG. Al. Chip design and false colored optical image of a superconducting qubit (inset). The chip comprises four supercon¬ 
ducting qubits Ql-4 (orange) made of aluminium and four niobium coplanar waveguide resonators Rl-4 (deep blue) coupled to 
input and output ports (red). The qubits have individual microwave drive lines (green) and flux bias lines (blue). 


ity band-pass filter (BP) and another circulator (C5) 
at the still stage were used to isolate the sample and 
JPD from higher-temperature noise. The transmitted 
signal is further amplified by a high electron mobility 
transistor (HEMT) at the 4.2 K stage and a chain of 
ultra-low-noise (ULN) and low-noise (LN) amplifiers at 
room temperature as shown in Fig. A2d. The amplified 
readout pulse is down-converted to an intermediate fre¬ 
quency (IF) of 25 MHz using an IQ mixer (Fig. A2e) 
and digitally processed by field-programmable gate ar¬ 
ray (FPGA) logic for real-time data analysis. 

Appendix B: Implementation of the XY gate 

The interaction between two qubits with degener¬ 
ate transition frequencies dispersively coupled to the 
same CPW resonator is described by the exchange cou¬ 
pling [33] J((ji<j 2 + 02 ”) which can also be written 

in terms of Pauli operators as ^ (erf erf + cr^erf). We acti¬ 
vate this interaction by tuning the transition frequency of 
qubit Q1 into resonance with qubit Q2 with a flux pulse 
(Fig. A3) for an interaction time r which we varied from 
0 to 60 ns. At the frequency of qubit Q2, we obtain a 


coupling strength J = —40.4 MHz from a fit to the spec¬ 
troscopically measured avoided crossing. To compensate 
overshoots of the flux pulse due to the limited bandwidth 
of the flux line channel, we use an inverted linear filter 
based on room-temperature response measurements of 
the flux line channel and in-situ Ramsey measurements 
of the residual detuning of qubit Q1 in the time interval 
from 0 to 2 /is after the flux pulse. 

Since the outcome of the XY gate depends strongly 
on the relative phase of the two-qubit input state, we 
have used the same LO signal for the upconversion of 
the single-qubit pulses acting on both qubits Q1 and Q2 
(green lines in Fig. A2a). Then the initial relative phase 
between the qubits is defined solely by the pulse sequence 
generated by the AWG and the cable lengths. In addi¬ 
tion, we choose the shape of the flux pulse that realizes 
the XY gate such that the dynamic phase acquired by 
qubit Q1 during the idle time and the rising edge of the 
flux pulse cancels any unwanted relative phase offset of 
the initial state. We satisfy this condition by tuning the 
frequency of Q1 to an intermediate level (buffer) for a 
fixed time of 16 ns before and after the XY gate (Fig. A3). 
A suitable buffer level is found by performing Ramsey- 
type experiments with a single XY gate while sweeping 
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FIG. A2. (a) Schematic of the experimental setup with complete wiring of electronic components inside and outside of the 
dilution refrigerator with the same color code as in Fig. Al. (b) Up conversion circuit for generating controlled microwave 
pulses, (c) Quantum limited parametric amplifier circuit to amplify readout pulses at base temperature, (d) Amplifiers used 
at room temperature just before down conversion of the signal, (e) Down conversion circuit (See text for details). 


the buffer amplitudes. This calibration procedure is car¬ 
ried out for each interaction length of the XY gate. The 
second buffer at the falling edge of the flux pulse is used 
to ascertain that the relative phase between the qubits 
after tuning qubit Q1 back to its original position is the 
same as the initial relative phase. 


Appendix C: Pulse scheme 

The quantum protocols for the digital quantum simu¬ 
lation of Heisenberg (Fig. A4a) and Ising spin (Fig. A4b) 


models were realized by sequences of microwave and flux 
pulses applied on qubit Q1 (red curves in Fig. A4) and 
qubit Q2 (blue curves in Fig. A4). The single-qubit rota¬ 
tions were implemented by 24 ns long Gaussian-shaped 
resonant DRAG [34, 35] microwave pulses and the XY 
gates were implemented using fast flux pulses. To avoid 
the effect of residual transient response of the flux pulse 
we have added a 40 ns + At waiting time after each flux 
pulse, with At being an adjustable idle time. We have 
chosen At such that the time difference between two ap¬ 
plications of the XY interaction is commensurate with 
the relative phase oscillation period of 5 ns, equal to the 
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7.14 GHz 
5.440 GHz 


5.240 GHz 


FIG. A3. Implementation of the XY gate. The transition fre¬ 
quency of qubit Q1 (red) is tuned into resonance with qubit 
Q2 (blue) for an interaction time r using a fast flux pulse. 
Before and after the flux pulse, a 16 ns long buffer is added 
at an intermediate level to cancel the dynamic phase accu¬ 
mulated by qubit Q1 relative to Q2 (grey area) during the 
evolution (see text). 


inverse frequency detuning 1/200 MHz. With these mea¬ 
sures we ensure that the gate can be used in a modular 
fashion, i.e. that a single calibration of the gate suffices 
for all gate realizations within the algorithm. The single¬ 
qubit phase gates were implemented by detuning the idle 
frequencies of each qubit with a square flux pulse. In 
the idle state, we observe a state-dependent qubit transi¬ 
tion frequency shift of 940 kHz due to the residual erf erf 
interaction. To decouple this undesired effect we have 
used a standard refocusing technique [36] implemented 
by two consecutive 7r pulses on qubit Q2 (magenta boxes 
in Fig. A4). In the end of each pulse sequence we per¬ 
form dispersive joint two-qubit state-tomography [37] by 
single-qubit basis transformations followed by a pulsed 
microwave transmission measurement through resonator 
Rl. 


SWAP gate [43, 44]. For a phase angle 7r, the Heisenberg 
interaction is an identity gate (Fig. A6b) with a process 
fidelity of 83.6 %. 


Appendix E: Error contributions 

The single-qubit gate fidelities measured by random¬ 
ized benchmarking [45-47] amount to 99.7%. The dom¬ 
inant contribution to the loss in fidelity originates from 
the two-qubit XY gates for which a process fidelity 
Jp^XY = 95.7% is obtained from process tomography 
averaging over all quantum phase angles. This indi¬ 
cates that the errors in the implementation of the XY 
gate limit the fidelity of the final state of the quantum 
simulation. To confirm this, we calculate the expected 
process fidelity for the Heisenberg and Ising protocol 
from the observed XY gate fidelity by assuming inde¬ 
pendent gate errors in all three steps. For the Heisen¬ 
berg (XYZ) model simulation neglecting the small single¬ 
qubit gate errors, we expect a mean process fidelity 
J P) xyz ~ 1 — 3(1 — T^xy) = 87.1 %, which is close to the 
observed value of 86.3%. For the Ising model simulation 
we expect a process fidelity of 7%,ising ~ 1—2n(l—T-^ xy)* 
From the relation P s = (dp p + l)/(d + 1) between state 
(P s ) and process fidelity (7^), we obtain the expected 
mean state fidelities of {93.1,86.2,79.4,72.5,65.6}% for 
n = 1 to 5 Trotter steps which compare well to the mea¬ 
sured state fidelites {91.7, 88.3, 82.2, 73.0, 60.7} %. 

To estimate the dominant source of systematic er¬ 
rors, we consider a model which includes relaxation (Ti) 
and dephasing (X 2 ) and state-dependent phase errors 
described by an effective erf erf term with interaction 
strength J z . In addition, we include an extra offset in 
the single-qubit phase gate acting on qubit Q2 from cross 
talk of the flux pulses acting on qubit Q1 in each Trotter 
step. By fitting the final state predicted by this model 
to the observed states, we estimate an unwanted inter¬ 
action angle J z r z of approximately 2.3° and a constant 
phase offset of 4.6°. 


Appendix D: Process tomography 

We perform standard two-qubit process tomo¬ 
graphy [38, 39] of the XY gate and of the simulated 
isotropic Heisenberg (XYZ) model for a varying inter¬ 
action time r. Fig. A5 shows the process x matrices 
characterizing the XY gate for a quantum phase an¬ 
gle 7r/2 (Fig. A5a) and 7 r (Fig. A5b) corresponding to 
a yjiSWAP gate [29, 40] and an iSWAP gate [41, 42] 
with process fidelities of 97.8% and 95.3%, respectively. 
Heisenberg interaction with a quantum phase angle tt/2 
leads to a SWAP gate (Fig. A6a) with a process fidelity of 
86.1%. While the SWAP gate belongs to the two-qubit 
Clifford group, there is no natural interaction in stan¬ 
dard circuit QED architecture to directly implement the 























FIG. A4. Pulse sequences that are applied on qubit Q1 (red) and qubit Q2 (blue) to implement the Heisenberg (a) and Ising 
spin (b) models. The Gaussian-shaped DRAG microwave pulses are applied to the charge lines of the respective qubits to 
implement single-qubit rotations R%, y about the x or y axis of the Bloch vector by an angle 0. Each sequence starts with the 
preparation of an initial state (green boxes) and ends with microwave pulses for basis rotations to perform state-tomography 
(yellow boxes). The microwave pulses marked with magenta boxes are used for refocussing. The black vertical bars with the 
two dots in panel (b) indicate that the enclosed pulse sequence is repeated n times. The XY gates are realized by applying flux 
pulses to the flux line of qubit Q1 for a time r/n. The phase gates Rt^ n are implemented by detuning the transition frequency 
of each qubit from their idle frequencies applying flux pulses for a time r/n. The numbers stated below the pulses on qubit Q1 
represent timescales in ns. 
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(a), XY (2 |J|t = tt/2) 
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(b), XY (2|J|x = 7i) 


95.3 % 
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FIG. A5. (a) Measured real and imaginary part of the XY process x matrix (Re x 5 Im x)> i n the basis {/ = identity , X — 
cr x , Y = —i(j y , Z — a z }, describing the mapping from any initial state to the final state for a quantum phase angle of 
2\J\r = 7t/2. The dashed wire frames represent the theoretically optimal matrix elements and the colored bars represent 
measured positive (blue) and negative (red) matrix elements. The fidelity of the experimentally observed process with respect 
to the ideal process is indicated in the black boxes, (b) As in (a) for a phase angle 7 r. 
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FIG. A6. (a) Measured real and imaginary part of the Heisenberg (XYZ) process x matrix (Re x? I m x)> m the basis 

{/ = identity , X — cr x ,Y = —ia y , Z = a z }, describing the mapping from any initial state to the final state for a quantum 
phase angle of 2| J\r m 7r/2. The dashed wire frames represent the theoretically optimal matrix elements and the colored bars 
represent measured positive (blue) and negative (red) matrix elements. The fidelity of the experimentally observed process 
with respect to the ideal process is indicated in the black boxes, (b) As in (a) for a phase angle i r. 
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